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The multifractal (MF) distribution of the electrostatic potential near any conformally invariant 
fractal boundary, like a critical O(N) loop or a Q -state Potts cluster, is solved in two dimensions. 

The dimension f(6) of the boundary set with local wedge angle 8 is f(6) = j — ^jSp e%Z-e) > with 
c the central charge of the model. As a corollary, the dimensions -Dep of the external perimeter 
and .Dh of the hull of a Potts cluster obey the duality equation (Dep — l)(Dn — 1) = 4. A related 
covariant MF spectrum is obtained for self-avoiding walks anchored at cluster boundaries. 

PACS numbers: 02.30.Em, 05.45.Df, 05.50.+q, 41.20.Cv 



Fractals are by now a well-recognized and studied sub- 
ject, yet they still pose great challenges in mathemati- 
cal physics: to establish a complete theory of diffusion 
limited aggregates (DLA), or a rigorous theory of the 
tenuous fractal structures arising in critical phenomena, 
to name only two. Classical potential theory, i.e., that 
of the electrostatic or diffusion field near such random 
fractal boundaries, whose self-similarity is reflected in a 
multifractal (MF) behavior of the potential, is almost a 
Terra Incognita. In DLA, the potential, also called har- 
monic measure, actually determines the growth process 
and its scaling properties are intimately related to those 
of the cluster itself JjJ. In statistical fractals, the Lapla- 
cian field is created by the random boundary, and should 
be derivable, in a probabilistic sense, from the knowledge 
of the latter. A singular example was studied in Ref. 
0j , where the fractal boundary, the "absorber" , was cho- 
sen to be a simple random walk (RW), or a self-avoiding 
walk (SAW), accessible to renormalization group meth- 
ods near four dimensions. 

In two dimensions (2D), conformal field theory (CFT) 
has lent strong support to the conjecture that statistical 
systems at their critical point produce conformally in- 
variant (CI) fractal structures, examples of which are the 
continuum scaling limits of RW's, SAW's, critical Ising 
or Potts clusters, etc. Given the beautiful simplicity of 
the classical method of conformal transforms to solve 2D 
electrostatics of Euclidean domains, perhaps a universal 
solution is possible for the planar potential near a CI 
fractal. 

A first exact example has been recently solved in 2D 
for the whole universality class of random or self-avoiding 
walks, and percolation clusters, which all possess the 
same harmonic MF spectrum || (see also B). In this 
Letter, I propose the general solution for the potential 
distribution near any conformal fractal in 2D. The exact 
multifractal spectra describing the singularities of the po- 
tential, or, equivalently, the distribution of wedge angles 
along the boundary, are obtained, and shown to depend 



only on the so-called central charge c, a parameter which 
labels the universality class of the underlying CFT. I de- 
vise conformal tools (linked to quantum gravity), which 
allow the mathematical description of random walks in- 
teracting with CI fractal structures, thereby yielding a 
complete, albeit probabilistic, description of the poten- 
tial. The results are applied directly to well-recognized 
universal fractals, like O(N) loops or Potts clusters. In 
particular, a subtle geometrical structure is observed in 
Potts clusters, where the external perimeter (EP), which 
bears the electrostatic charge, differs from the cluster's 
hull. Its fractal dimension Dep is obtained exactly, gen- 
eralizing the recently elucidated case of percolation jj| . 

Harmonic Measure and Potential. Consider a single 
(conformally invariant) critical random cluster, gener- 
ically called C. Let H (z) be the potential at exte- 
rior point z £ C, with Dirichlet boundary conditions 
H (w S dC) = on the outer (simply connected) bound- 
ary dC of C, and H(w) — 1 on a circle "at 00", i.e., of a 
large radius scaling like the average size R of C. From a 
well-known theorem due to Kakutani Q, H (z) is iden- 
tical to the harmonic measure, i.e, the probability that 
a random walker launched from z, escapes to 00 without 
having hit C. The multifractal formalism J7|— |lC|] charac- 
terizes subsets dC a of boundary sites by a Holder expo- 
nent a, and a Hausdorff dimension /(a) = dim(9C Q ), 
such that their potential locally scales as 



H(z^we dC a ) « (|z - w\/R) a , 



(1) 



in the scaling limit a -C r = \z — w\ « JJ, with a the 
underlying lattice constant. In 2D the complex potential 
ip(z) (such that the electrostatic potential H(z) = ^tcp(z) 
and field |E(z)| = |<//(z)|) reads for a wedge of angle 0, 
centered at w: (p(z) = (z — w)*/ e . By Eq. ([!]) a Holder 
exponent a thus defines a local angle = ir/a, and the 
(purely geometrical) MF dimension f(0) of the boundary 
subset with such 9 is f{9) = f(a = tt/9). 

Of special interest are the moments of H, averaged 
over all realizations of C 
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where 9C(r) is shifted a distance r outwards from dC, 
and where n can be a real number. In the scaling limit, 
one expects these moments to scale as 



Z n « (r/i?) r 



(n) 



(3) 



where the multifractal scaling exponents r (n) encode 
generalized dimensions, D (n) — t (n) / (n — 1), which 
vary in a non-linear way with n f^-[lOt; they obey the 
symmetric Legendre transform r (n) + f (a) = an, with 
n = f (a) , a = r' (n). From Gauss' theorem Q r(l) = 
0. Because of the ensemble average (^), values of / (a) can 
become negative for some domains of a Q . This Letter 
is organized as follows: I first present the main results in 
a universal way, then proceed with their derivation from 
conformal field theory, and finally specify them for the 
O(N) and Potts models. 

Exact Multifractal Dimensions and Spectra. Each con- 
formally invariant random system is labelled by its cen- 
tral charge c, c < 1. The multifractal dimensions of a 
simply connected CI boundary then read explicitly: 
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The Legendre transform reads 
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Notice that 
any c, r'(n 

/(« = 1) = 
any simply 
also remark 
3-2L>(n) = 
bounds: Vc. 
< 9 < 2tt 
distribution 



the generalized dimensions D(n) satisfy, for 
= 1) = D(n = 1) = 1, or equivalently 
= 1, i.e., Makarov's theorem [ jn] , valid for 
connected boundary curve. From (|]J^) we 
a fundamental relation, independent of c: 
l/a= 9/tt. We also have the superuniversal 
D(oo) < D(n) < D(n*) = §, hence 
We arrive at the geometrical multifractal 
of wedges 9 along the boundary: 



Vn,i 



/» = /© = J 



25 



{■K-9Y 



12 6>(2tt 



(7) 



Remarkably enough, the second term also describes the 
contribution by a wedge to the density of electromag- 
netic modes in a cavity p"2| . The maximum of /(a) cor- 
responds to n — 0, and gives the dimension Dep of the 
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FIG. 1. Universal harmonic multifractal spectra f(a). The 
curves are indexed by the central charge c, and correspond 
respectively to: 2D spanning trees (c = —2); self-avoiding or 
random walks, and percolation (c = 0); Ising clusters or Q — 2 
Potts clusters (c = |); N = 2 loops, or Q = 4 Potts clusters 
(c = 1). The maximal dimensions are those of the accessible 
frontiers. The left branches of the various f(ct) curves are 
largely indistinguishable, while their right branches split for 
large a, corresponding to negative values of n. 



support of the measure, i.e., the accessible or external 
perimeter as sup Q /(a) = /(a(0)) = D(0): 
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This corresponds to a typical wedge angle 9 
7r/a(0) = 7r(3- 2£>ep)- 

The multifractal functions / (a) — a = f(6) — ? are 
invariant when taken at primed variables such that 2-7T = 
(9 + 9' = ^ + ^7 , which corresponds to the complemen- 
tary domain of the wedge 9. This condition reads also 
D(n) + D(n') = 2. This basic symmetry, first observed 
and studied in |l3| for the c = result of || , is valid for 
any conformally invariant boundary. 

In Fig. 1 are displayed the multifractal functions /, 
Eq. (^), corresponding to various values of —2 < c < 1, 
or, equivalently, to a number of components N G [0,2], 
and Q e [0,4] in the O(N) or Potts models (see below). 
The generalized dimensions D(n) (not shown) appear 
as quite similar for various values of c, and a numeri- 
cal simulation would hardly distinguish the different uni- 
versality classes, while the /(a) functions do distinguish 
these classes, especially for negative n. The singularity 
at a = |, or 6 = 2ir, in the multifractal functions /, 
or /, corresponds to boundary points with a needle lo- 
cal geometry, and Beurling's theorem p4| indeed insures 
the Holder exponents a to be bounded below by | . This 
corresponds to large values of n, where, asymtotically, for 
any universality class, Vc, linin^oo D(n) = h. 

The right branch of / (a) has a linear asymptote 
lim Q ._ > .oo / (a) /a — n* = —(1 — c)/24. The limit mul- 
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tifractal spectrum is obtained for c = 1, which exhibits 
an exact example of a left-sided MF spectrum, with an 
asymptote / (a— > oo, c = 1) — > § (Fig. 1). It corre- 
sponds to singular boundaries where / (9 — > 0, c = 1) = 
| = -DePi where the external perimeter is dominated 
by " fjords" , with typical angle = 0. 

The a — > oo behavior corresponds to moments of low- 
est order n — > n* , where D(n) reaches its maximal value: 
Vc, D(n*) = |, common to all simply connected, confor- 
mally invariant, boundaries. This describes almost inac- 
cessible sites: define Af (H) as the number of boundary 
sites having a given probability H to be hit by a RW 
starting at infinity; the MF formalism yields, for H — > 0, 
a power law behavior Af (if) \h^o ~ H~^ 1+n > with an 
exponent 1 + n* = 2|±£ < 1. 

Conformal Invariance and Quantum Gravity. Let me 
now give the main lines of the derivation of exponents 
7"(n), hence f(a), by generalized conformal invariance. 
By definition of the ii-measure, n independent RW's, or 
Brownian paths B in the scaling limit, starting at the 
same point a distance r away from the cluster's hull dC, 
and diffusing without hitting dC, give a geometric repre- 
sentation of the n th moment, H n , in Eq.(|^) for n integer. 
Convexity yields analytic continuation for arbitrary n's. 
Let us introduce the notation AAB for two random sets 
conditioned to traverse, without mutual intersection, the 
annulus T> (r, R) from the inner boundary circle of radius 
r to the outer one at distance R, and A V B for two in- 
dependent, thus possibly intersecting, sets |j| . With this 
notation, the probability that the Brownian paths and 
cluster are in a configuration dC A (V£>)™ = dC A n, is 
expected to scale for R/r oo as 



V R (dC An)* (r/R) 



x(n) 



(9) 



where the scaling exponent x (n) depends on n. In terms 
of definition (^) , the harmonic measure moments (|^) sim- 
ply scale as Z n w R 2 Vr (dC A n) (||| j which, combined 
with Eq. (||), leads to t (n) = x (n) — 2. 

To calculate these exponents, I use a fundamental map- 
ping of the conformal field theory in the plane R 2 , de- 
scribing a critical statistical system, to the CFT on a fluc- 
tuating abstract random Riemann surface, i.e., in pres- 
ence of quantum gravity Jl5| , ^6[ . Two universal functions 
U, and V, depending only on the central charge c of the 
CFT, describe this map: 



U(x) =x x —L, V{x) = \^—L 
1 — 7 4 1 — 7 



(10) 



with V (x) = U (| (x + 7)) [0. The parameter 7 is the 
string susceptibility exponent of the random 2D surface 
(of genus zero), bearing the CFT of central charge c Jl5| ; 
7 is the solution of c = 1 — 67 2 (1 — 7) _1 ,7 < 0. 

For two arbitrary random sets A, B, with boundary 
scaling exponents in the half-plane x (A) , x (B) , the scal- 



ing exponent x (A A B) , as in 
ture §l7l 



has the universal struc- 



x(AAB) =2V [if- 1 {x (A)) + U- 1 (x (B))] , (11) 
where U~ x (x) is the positive inverse function of U 



ir 



(x) = - (V4(l-7)z + 7 2 +7) 



(12) 



C/ _1 (x) is, on the random Riemann surface, the bound- 
ary scaling dimension corresponding to x in the half- 
plane R x R + , and the sum of Z7 -1 functions in Eq. ( |Tl| ) 
is a linear representation of the product of two "bound- 
ary operators" on the random surface, as the condition 
AAB for two mutually avoiding sets is purely topological 
there. The sum is mapped back by the function V into 
the scaling dimensions in R 2 [[l7j . 

For the harmonic exponents x(n) = x (dC A n) in 
we use (pT[). The external boundary exponent x (dC) 
obeys U (x) = 1 — 7, which I derive either directly, 
or from Makarov's theorem: ^(n = 1) = 1 ||l7||. The 
bunch of n independent Brownian paths have simply 
x ((V£) n ) = n, since x(B) = l [§. Thus I obtain 



x{n) = 2V{l -"t + lJ- 1 (n)) . 
This finally gives from ©(jl2[) r(n) = x(n) - 2 



(13) 



T ( n ) = ^ n 



1) 



12-7 
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[^4(1 - 7 )n + 7 2 - (2 - 7)] 



from which Eq. (^) follows, QED. 

This formalism immediately allows generalizations. 
For instance, in place of n random walks, one can con- 
sider a set of n independent self-avoiding walks V , which 
avoid the cluster fractal boundary, except for their com- 
mon anchoring point. The associated multifractal expo- 
nents x (dC A (VP)™) are given by (^3|), with the argu- 
ment n in l/- 1 simply replaced by x ((VP)") = nx (V) = 
In These exponents govern the universal multifractal 
behavior of the moments of the probability that a SAW 
escapes from C. One then gets a spectrum / such that 
f(a = x (V) tt/60 = / (a = w/0) = f(0), thus unveiling 
the same invariant underlying wedge distribution as the 
harmonic measure, QED. 

0{N) and Potts Cluster Frontiers. The 0(N) model 
partition function is that of a gas Q of self- and mutu- 
ally avoiding loops on a given lattice, e.g., the hexag- 
onal lattice Hi: Z 0{N) = J2g K Ub N Mp , with K and 
N two fugacities, associated respectively with the total 
number of occupied bonds Mb , and with the total num- 
ber Afp of loops, i.e., polygons drawn on the lattice. For 
N € 0.[— 2, 2], this model possesses a critical point (CP), 
K c , while the whole "low-temperature" (low-T) phase, 
i.e., K c < K, has critical universal properties, where the 
loops are denser that those at the critical point |l8| . 

The partition function of the Q-state Potts model on, 
e.g., the square lattice, with a second order critical point 
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for Q £ [0,4], has a Fortuin-Kasteleyn representation at 
the CP: Zpotts = J2u(C) Q^ F where the configurations 
U(C) are those of reunions of clusters on the square lat- 
tice, with a total number Mp of polygons encircling all 
clusters, and filling the medial square lattice of the origi- 
nal lattice (l^Jl^l ■ Thus the critical Potts model becomes 
a dense loop model, with a loop fugacity N = Q%, while 
one can show that its tricritical point with site dilution 
corresponds to the O(N) CP §0|. The 0(N) and Potts 
models thus possess the same "Coulomb gas" represen- 
tations N = s/Q = -2 cos 71-5, with g £ [1, §] for 
the 0{N) CP, and g € [§, 1] for the low-T O(N), or criti- 
cal Potts, models; the coupling constant g of the Coulomb 
gas yields also the central charge: c = 1 — 6(1 — g) 2 1 9- 
The above representation for N = \[Q £ [0,2] gives a 
range of values — 2 < c < 1; our results also apply for 
c£ (—oo, —2], corresponding, e.g., to the O (N £ [—2,0]) 
branch, with a low-T phase for g £ [0, |], and the CP for 

<?e[f,2]. 

The fractal dimension Dep of the accessible perimeter, 
Eq. d§|), is, like c(g) = c(g^ 1 ), a symmetric function 

^ep = i + - g- 1 ) + \g$(i - g), (14) 

where $ is the Heaviside distribution, thus given by two 
different analytic expressions on either side of the sepa- 
ratrix g = 1. The dimension of the hull, i.e., the com- 
plete set of outer boundary sites of a cluster, has been 
determined for O(N) and Potts clusters pl| j, and reads 
£>h = 1 + iff 1 ) f° r the entire range of the coupling 
constant g £ [|,2]. Comparing to Eq. (jljj), we there- 
fore see that the accessible perimeter and hull dimen- 
sions coincide for g > 1, i.e., at the O(N) CP (or for 
tricritical Potts clusters), whereas they differ, namely 
T>ep < Djj, for g < 1, i.e., in the O(N) low-T phase, 
or for critical Potts clusters. This is the generalization to 
any Potts model of the effect originally found in perco- 
lation |^2[ . This can be directly understood in terms 
of the singly connecting sites (or bonds) where fjords 
close in the scaling limit. Their dimension is given by 
D sc = 1 + \g- x - \g §|. Thus, for critical O(N) 
loops, g £ (1,2] and Z?sc < 0, so there exist no closing 
fjords, which explains the identity: -Dep = -Dh! whereas 
Dgc > 0,3 £ [hi 1) for critical Potts clusters, or in the 
O(N) low-T phase, where pinching points of positive di- 
mension appear in the scaling limit, so that -Dep < £>h 
(Table 1). I then find from Eq. @, with g < 1: 

(D EP - 1) (T> H - 1) = \. (15) 

The symmetry point Dep = -Dh = § corresponds to g = 
1, N = 2, or Q = 4, where, as expected, the dimension 
T*sc = of the pinching points vanishes. 

For percolation, described either by Q = 1, or by the 
low-T 0(N = 1) model, with g — I, we recover the result 
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TABLE I. Dimensions for the critical Q-state Potts model; 
Q = 0, 1, 2 correspond respectively to spanning trees, perco- 
lation and Ising clusters. 



Dep = |, recently derived in M. For the Ising model, 
described either by Q = 2, g = j, or by the 0(N = 1) 
CP, g' = g^ 1 — |, we observe that the EP dimension 
Dep = ^ coincides, as expected, with that of the critical 
0(N — 1) loops, which in fact appear as EP's. This is 
a particular case of a further duality relation between 
the critical Potts and CP 0(N) models: T> E p (Q(g)) = 
D K (0(N(g')))Jox g' =g-\g<l. 
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